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ABSTRACT 

Fast Quasi-Periodic Oscillations (QPOs, frequencies of ~ 20-1840 Hz) have been 
recently discovered in the ringing tail of giant flares from Soft Gamma Repeaters 
(SGRs), when the luminosity was of order lO'*' - 10"^^^ erg/s. These oscillations per- 
sisted for many tens of seconds, remained coherent for up to hundreds of cycles and 
were observed over a wide range of rotational phases of the neutron stars believed 
to host SGRs. Therefore these QPOs must have originated from a compact, virtually 
non-expanding region inside the star's magnetosphere, emitting with a very moder- 
ate degree of beaming (if at all). The fastest QPOs imply a luminosity variation of 
AL/ A? ~ 6 X lO'^^ erg s"^, the largest luminosity variation ever observed from a com- 
pact source. It exceeds by over an order of magnitude the usual Cavallo-Fabian-Rees 
(CFR) luminosity variability limit for a matter-to-radiation conversion efficiency of 
100%. We show that such an extreme variability can be reconciled with the CFR limit 
if the emitting region is immersed in a magnetic field > 10'^ G at the star surface, 
providing independent evidence for the superstrong magnetic fields of magnetars. 



Subject headings: stars: magnetic fields — stars: neutron — stars: individual(SGR 
1806-20) 
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1. Introduction 



Soft Gamma Repeaters are a small class of galactic sources of X and soft gamma radiation. 
They have spin periods of ~ 5 10 s, display a secular spin-down with timescales of ~ 10"^ 10^ yr 
and do not possess a companion. Unlike radio pulsars, the rotational energy loss of SGRs is a factor 
of 10 -h 100 too small to ex plain their persistent emission, typically ~ 10^^ 10^^ erg/s (see e. g. 
Woods & ThompsonI hoo4\ \ Like Anomalous X-ray Pulsar (AXPs, iMereghetti & Stellal dlQQSh . 
with whom they share a number of properties, SGRs are believed to host magnetars, neutron stars 
the emi ssion of which is powered by the decay of their superstrong (internal) magnetic field {B > 

G. buncan & ThompsonI (Il992h : [Thompson & Duncanl(ll993h . 



The name-defining characteristic of SGRs is that they show periods of activity in which re 
current short bursts are emitted, with peak luminosities of ~ 10^^ 10"^^ ergs s"' and sub-secom 
durations. The characteristics of the three giant flares observed so far in about 30 yr of mon 
itoring are much more extreme. Their initial, tenths-of-seconds-long spike releas es enormous 
amounts of energy, ~ 10"*"^ ergs in the 1979 March 5 eve nt from SGR 0526-66 (iMazets et al. 



19791) and the 1998 August 27 event from SGR 1900-1-14 (,Hurlev et all 1 19991 : Feroci et al.l[ l999l) 
and as much as 5 x 10^^ erg in the 2004 December 27 event from SGRl 806-20 d Hurlev et al l EoOS: 
Mereghetti et allboosljpalmer et allboOS^Terasawa et allboosl) . After the initial spike, giant flares 
display a very bright ringing tail lasting hundreds of seconds and releasing a total energy of about 
~ 10"^"^ ergs. The emitted spectrum i s roughly thermal, w ith a blackbody equivalent temperature of 
~ 5 keV in the case of SGRl 806-20 iHurlev et all toosl) . 



The highly super-Eddington luminosities of the recurrent bursts of SGRs and especially of 
their giant flares make models involving accretion energy not viable. According to the magne- 
tar model, the emission of SGRs (and A XPs) draws from their extremely high magnetic fields 
(|Thompson & Duncanlll995L 1 1 99611200 ih . Within this model the neutron star interior is character- 
ized by a wound-up, mainly toroidal magnetic field configuration with Bt > 10'^ G. A less intense 
(mainly poloidal) field emerges o ut of the star magnetosphere, cau sing spin-down via rotating 
dipole losses at the observed rate ([Thompson & Duncanlll993L 12001 ) (dipole B-field strength of 
Bel ~ 7.8 X 10^"^ G). Energy propagates to the neutron star magnetosphere through Alfven waves 
driven by local "crust-quakes" and giving rise to recurrent bursts with a large range of amplitudes. 
Large-scale rearrangements of the internal magnetic field or catastrophic instabilities in the mag- 
netosphere are invoked to explain the sudden release of very large amounts of energy that occurs 
in giant flares ([Thompson & Duncanll200ll : lLyutikovll2003|) . A fireball of plasma expanding at 
relativistic speeds breaks out of the star's magnetosphere, causing the initial sub-second spike of 
giant flares. The ringing tail that follows results most likely from the part of the fireball that re- 
mains trapped in the star's magnetosphere. The energy release in the ringing tail yields a limit for 
the external field of magnetars ( ^ 10'"^ G) in agreement with the values inferred from spin-down 
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dipole losses (|Thompson & Duncanlll995Ll2001|) . wh ile an analysis of the i nitial time scales gives 
evidence in favor of the crustal cracking mechanism ([Schwartz et al. Il2005|) . 



To confirm this model it is thus essential to me asure the surface m agnetic field, which com- 
plements the measurement of the dipole component (IWoods et al.l (l2002h ). to include higher order 
multipoles. In this paper we present an interpretation of the QPOs observed in SGR 1806-20 which 
provides a lower limit (B ^ B^pt 4.4 x 10^^ G) to the surface field. The argument is based upon 
a seldom-used constraint ( dCavallo & Reeslll978l : iFabian & Reeslll979l : lFabianlll979l) ) that puts a 
very strong upper limit on the time-scale on which significant luminosity variations can take place. 
In the next section we discuss the relevant observations; in Section 3 we will re-derive the limit 
and show that it is largely violated by QPOs in SGR 1806-20; after vainly trying to circumvent it, 
we will show that it can be reconciled with observations only if B Bq. 



2. Quasi Periodic Oscillations in Giant Flares of SGRs 



Recent studies led to the discovery that the X-ray flux of th e ringing tail of S GRs' giant 
flares is characterized by fast Quasi Periodic Oscillations, QPOs (^ Israel et aL 2005 ). Different 
QPO modes were detected, some of which were excited simultaneously. The ringing tail of the 
December 2004 e vent from SGR 1806-20 disp layed clear QPO signals at about 18, 30, 93, 150, 
625 and 1840 Hz Jstrohmayer & WattsI ho06^ ). Similarly, QPOs around frequencies of 28, 54, 
84 and 155 Hz were detec ted during the ringing tail of the 1998 giant flare of SGR 1900-1-14 
(|Strohmayer & WattsI l2005h. while hints fo r a signal at ~43 Hz were found in the March 1979 
event from SGR 0526-66 (|Barat et al.ll 19831) . These QPOs show large variations of the amplitude 
with time and, especially, of the phase of the spin modulation in the giant flare's tail. 

The similarity in some of the QPO modes and frequencies across different SGRs suggests that 
the production mechanism is the same. A likely interpretation involves the excitation of neutron 
s tar oscillation r nodes, whose expected eigenf requencies match some of the observed QPOs peaks 
( (|Duncanlll998l ; llsrael et al.ll2005l ; |Pirdl2005h ). Infact, if giant flares result from large scale frac- 
turing of the crust induced by instabilities of the internal magnetic fields, then the excitation of 
crustal and (possibly) global neutron star modes is to be expected ( (ILevinll2006h ). Regardless of 
the exact mechanism driving these oscillations we are concerned here with the extremely large and 
fast luminosity variations of the QPOs. 

We concentrate on the signals with the largest luminosity tim e derivative, i.e. the 625 and 
1840 Hz QPOs from SGR 1806-20 dstrohmaver & WattsI J2006a b). The power spectrum peaks 
through which these QPOs are revealed, are a few Hz wide, testifying that their signal remained 
coherent for hundreds of cycles. The signal shape must be close to sinusoidal, as evidenced by 



-4- 



the absence of detectable harmonic signals. The 625 and 1840 Hz QPOs were detected only in 
a ~ 50 s long interval of the ringing tail, about 200 s after the initial spike, and were especially 
prominent over a ~ 140 deg interval in rotational phase. The QPO amplitude reached a maximum 
in this phase interval over two consecutive rotation cycles: for both signals the rms amplitude 
was arms ~ 18%. Approximating the QPOs with sinusoids, we estimate their highest luminosity 
derivative as AL/Af = l^^^irLarmsi^QPo, with uqpq the QPO frequency. Here L ~ lO'^' ergs s"^ is 
the luminosity in the relevant section of the ringing tail (for the likely source distance of 15 kpc). 
This gives AL/At ~ 1 x lO^-* and 3 x lO^^ ergs s'^ for the 625 and 1840 Hz QPOs, respectively. 
The effects of beaming might decrease these values somewhat, but not by a very large factor. In 
fact, these QPO signals were observed over a large interval of rotational phases (about ~ 140 deg), 
translating into approximately the same azimuthal range of emission angles. It is natural to assume 
a comparably large angular spread in latitude (unless the neutron star rotation axis is very close to 
our line of sight, which is unlikely given the large amplitude of the rotation modulation and the 
size of the emission region, see below). This gives a solid angle of order ~ tt ster. Adopting this 
beaming factor the luminosity derivatives above reduce to ~ 2 x lO'^^ and 6 x 10"^^ ergs s"^; this are 
the values that we adopt in the following discussion. 

We stress here, because this is essential to our argument (to be presented shortly), that, to- 
gether with the QPOs, in the ringing tail a strong modulation at the star's spin period is clearly 
present, similar in relative amplitude and shape to the modulation observed when the source is in 
its quiescent state: this (together with the lack of significant amounts of beaming) indicates that 
the emission in the ringing tail originates from a region that remains stably anchored to the star's 
magnetosphere, and thus that relativistic bulk motions are not present at this stage of the flare. 

The blackbody temperature and luminosity in the ringing tail translates into a lower limit on 
the size of the emitting region of about ~ 30 km, i.e. substantially larger than the neutron star. On 
the other hand the black body-like spectral shape testifies that the emitting region is optically thick 
(or at least effectively thick), implying a scattering optical depth ^ 1. We remark that the size 
estimate, ~ 30 km, will play an important role in the following. 



3. The Cavallo-Fabian-Rees Variability Limit 



There is a well-known limit on the rate of c hange of the luminosity of any given source , 
which we briefly summa rize here dCavaUo & ReesI (Il978l) : Ipabian & ReesI (Il979l) : Ipabianl (Il979l) . 
see also Lawrence! (| 1 9 80|) : iHoshil (| 1 9 89|) ) . Suppose a source undergoes a large luminosity variation 
on a time-scale At, and there is a luminosity variation, over this time-scale, AL. The total energy 



-5- 



released within A? is related to the total mass within the source dimension R by 



An 

ALAt = 7] —R^nnipC^ , ( 1 ) 



where n is the average baryon density, and rj is the energy extraction efficiency. The time- scale At 
must obviously exceed the time over which the photons manage to diffuse out of the source: 

At>-{\+TT) (2) 

c 

where tt = ajuR is the Thomson optical depth and aj the Thomson cross-section. Eliminating R 
from the first two equations, we find 

3^AL(^ ,3) 

Regarded as a function of Tj, the above has a minimum for tj = 1, which yields the limit 

3 (Tr AL 

> . (4) 

The above limit is remarkable in that it is independent of both R and n, or any combination thereof: 
only the dependence on AL, a directly observable quantity, is left. The Cavallo-Fabian-Rees, CFR, 
limit thus writes 

AL/At <ri2x lO'^^erg s'^ . (5) 



The 625 and 1840 Hz QPO signals from SGRl 806-20 exceed the CFR limit by about an order 
of magnitudeQ: the largest value found in the previous section is AL/At = 6 x 10"^^ erg s~'^, which 
is a whole factor 30 /rj larger than the CFR's limit. 

In order to appreciate how hard it is to circumvent this limit, notice the in situ re- acceleration 
of electrons does not help, because it does not change the energies reached by protons; if protons 
were to escape, leaving electrons behind to be re-accelerated at will. Coulomb forces would quickly 
make the escape of protons impossible. Nor will having relativistic protons help, as one might 
think, since this would imply energies per electron to be radiated ~ jmpC^ ^ rfnipC^, because the 
energy transfer from protons to electrons is too slow. To show this, let us consider what happens 
when protons transfer promptly to the electrons their energy gain: protons may then be relativistic, 
in which case the maximum energy which can be extracted from each of them is 'jnipC^, rather than 



'We note that also the slower - 90-150 Hz QPOs signals in the giant flares' tail of SGR1806-20 and SGR1900+14, 
exceed the CFR limit, though by a smaller factor. 
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the more sedate rjmpC^. This however appears like an unlikely way out, because, even admitting 
that the electrons radiate very promptly their internal energy, the time-scale on which protons 
manage to transfer to the electrons their internal energy is much longer than the above limit. To 
see this, we idealize the situation as one where electrons are cold {i.e. Newtonian) as a result of 
their short cooling timescales, while the protons are still relativistic. The energy transfer rate is 
given by the usual formula 

dEp 2ne^n^ ( , ,2m^v^V7, 




ln(— f^) + l-2- , (6) 



where is the electron number density, Up = y^Ane^ne/nie the plasma frequency, v ~ c the proton 
speed, and W ~ 2E^/{meC^) is the maximum energy transfer in the not too extreme limit, 7p < 
mp/nie (but considering the opposite limit, 7 > mp/nie would change the argument very little since 
W only appears as an argument to a logarithm). 

The energy transfer timescale is of course tg = Ep/(c dEp/dx), which, for this whole idea to 
work, must be shorter than At. The condition < At can also be rewritten as, in the limit ~ c^: 



7 3aT L e n 



< — — — (7) 

ln( ) - 1 mpC^T] nipm.c^ 

h ujf, 

Inserting the numerical values for the luminosity of SGR 1806-20 during the ringing tail and for 
the typical magneto spheric density for a pulsar with a normal magnetic field we find 

^<10-8l . (8) 

^ ~ 3 X 1041 ergs-^ 10"' cmr^ ^ ' 

In order to reconcile the CFR limit with observations, we need 7 ^ 30, corresponding to a magne- 
tospheric density n ^ 3 x 10^^ cm~^ . While it is certainly true that plasma outside pulsars need not 
be charge-separated, still this density would exceed the minimum (i.e., charge separated) value by 
more than 9 orders of magnitude, making it unlikely that such a plasma may exist. 

This argument must be modified in the presence of a pair plasma, where pair creation pro- 
cesses can easily lead to a st rong increase in n„ ocT^ . Still, we know that astrophysical pair plas- 



mas are thermally regulated (|Pietrini & Kroliklll995h . so that their temperatures always lie around 



r < 1 MeV. At these temperatures, ^ 10^^ cm~^, but we also know that Ep <C nipC^, so that 
the limit of eq. |4] still holds, provided the total energy density is still dominated by the rest mass 
of the baryons. When instead the total energy density is dominated by the pairs (in other words, 
when little or no baryons are admixed), eq. |2]obviously still holds (with tj now indicating the total 
optical depth due to the pairs and possibly to photons as well), while eq. [Umust be rewritten as 

ALAt = rj^R^nm.c^- (9) 
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where we used the fact that 7 2 for the electrons. The new limit for the variability is 

A? > ^— = 134 s-—-— , (10) 

which is even more stringent than eq. IH 

The CFR limit might fail if coherent phenomena are involved, but we are unaware of any 
coherent mechanism working in the X-ray region of the spectrum where QPOs are observed, 
rather than the usual radio band. 

The classic way to circumvent the CFR limit is by means of relativistic effects. After all, 
GRBs' light-curves often display millisecond variability, and even when they don't they turn on on 
timescales of a second or less, reaching L ^ 10^" erg s~^: this clearly violates eq. |4]by eight to ten 
orders of magnitude (for i] = 1). This is due to a combination of relativistic aberration, blueshift 
and time contraction, which can drastically increase the luminosity and decrease the variability 
timescale in the observer's frame. Also BL Lac objects an d Quasars have lumi nosity derivatives 



that exceed by orders of magnitude the CFR limit (see e.g. (iBassani et al.lll983|) ). even after their 
luminosity is corrected for beaming. This explanation is probably correct also for the giant flare of 
SGR 1 806-20, where values as high as AL/ At ~ lO'*^ ergs s~'^, for isotropic ejectior§ are observed: 
the observation of a radio halo expanding at relativistic speed testifies that the observed variability 
in the initial spike of the giant flare is most likely magnified by relativistic effects as well. 

However, this explanation is not suitable for the QPOs in the ringing tail of the same flare: 
in fact, the presence of modulations at the spin period, with similar pulse shape and amplitude 
to those during quiescent periods, makes the existence of relativistic bulk motions very unlikely 
during this phase. Also, one should remember that the photosphere size is estimated to be 30 km 
at all times during the ringing tail, which implies that emission during this phase is due to material 
stably anchored to the pulsar. Lastly, the observed temperatures (T ^ 5 keV) testify to the gas 
having thermal speeds much below the escape velocity at that radius. For these three reasons, we 
deem relativistic effects an unlikely explanation for the violation of eq. HI 

The last way in which the CFR limit can fail is the one we propose: the scattering cross-section 
may differ from Thomson's because of the presence of a strong magnetic field B, exceeding the 
quantum value Bq = rn^c^ /(eh) = 4.4 x lO^^G. In this case, the scattering cross section for the 
ordinary (Q) and extrao r dinary (E) modes, and for the conversion of photons into the other state. 



are given by iMeszaros I (|l992h . when the dielectric tensor is dominated by vacuum polarization 
effects, as: 

3 

dao^o = -^ctt sin^ 6 sin^ 6'dcos 6' 



^Beaming in the ejection, which is presently unknown, could decrease this value substantially. 
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. 3 / 

8 \meC^B 

A 3 / 

doE^o = c 

8 \ MeC-^B 



cos^ 9d COS 9' 



cos^e'dcosO' 

2 



dcosO' 



(11) 



where e is the photon energy, and 9 and 6*' are the angles between the photon momentum before 
and after the diffusion, respectively, with the direction of the magnetic field. 

This equation shows immediately that, for photons emitted in the extraordinary mode, the 
cross section is reduced, with respect to the Thomson value, by a factor ^ (eB^/(meC^B))^. Thus, 
to bring the observed value, AL/ At = 6 x 10'^^ erg s~^, in agreement with eq. HI we just need to 
have {eBq/{meC^Byf < 77/30; here we take for e the value e ~ 14 keV , which is the peak of the 
Planck distribution for the observed temperature T = 5 keV. So our conclusion is that the QPOs' 
luminosity variation agrees with the CFR limit provided 



B>1.5BJ — 

V 



1/2 



6.6x lO^^G. 



(12) 



A technical comment is in order at this point: th e description of r adiation transfer in terms of 
separate modal propagation is not always adequate (ILai & Ho I (l2003h ). because of mode collapse. 
However this effect seems to be mostly relevant for even higher fields (5 > 7 x 10^^ G) than those 
derived here, which means that our naive treatment is probably justified. 

We now remark that this lower limit applies to the field close to, but not at the surface of SGR 
1806-20, because, as discussed in the previous section, emission from the ringing tail is generated 
within 30 km. At this distance from th e star surface, the dipole field inferred from pulsar spin-down 
(5 = 7.8 X 10^4 G. IWoods et aP (l2002l) ). is 2.9 x 10" G(R„J 10 km)\ which is smaller than the limit 
just derived, as expected if higher order multipoles are relevant in the star vicinity. 

Given the rapid decrease of the dipole field (and a fortiori of the other multipoles) with dis- 
tance from the star, the surface magnetic field must certainly satisfy 



B > 1.8 X 10" G 



10 km 



0.1 

V 



1/2 



(13) 



where a dipole-like radial dependence has been assumed between 30 km and the star surface, 
located at i?,,^; we shall refrain fror n making more elaborate hy pothesis about the structure of 
the magnetic field within R = 30 km ([Thompson & DuncanI (1200 ih ). because our aim is simply to 
provide a lower limit to the surface field, for which this minimum hypothesis (pure dipole) is fully 
adequate. 
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4. Discussion 



Our final result is the lower limit in eq. [131 This value is close to that of the dipole field, as 
inferred from pulsar spin-down, 5 = 7.8 x \0^^ G, but our estimate includes higher order multipoles, 
at least at 30 km. It is thus completely independent of the estimates from pulsar spin-down. We 
note that the limit on B derived by requiring that the magnetic field close to the star can prevent 
the escape of the "trapped fire ball" includes also the contribution of higher multipole components 



ipp 
in|( 



(|Thompson & DuncanI (120011) ). but gives a substantially lower magnetic field (B > 10 G) than 



the limit derived above. 

There is an easy way to test our interpretation of the failure of the CFR limit. The O- and 
E-mode photospheres of the trapped fireball are located at different heights in the star's magneto- 
sphere, as a result of the different electron scattering cross sections, which therefore sample regions 
of different B-field strengths and orientations. The polarization fraction and angle of the emerging 
X-ray flux should thus be modulated with the phase of the QPO signal as a result of the varying 
(relative) intensity of the O- and E-mode photon component. This in principle gives a clear test for 
the correctness of our interpretation, which might perhaps become verifiable in the future. 



Photons in the 0-mode suffer a strong Comptonization ([Thompson & DuncanI (|2001[) ). If 
the atmosphere were due to pure scattering, the ensuing photon distribution would differ from a 
blackbody at low photon frequency, oc instead of oc z/^, a result due to photon number con- 
servation. At first sight, one might think that there are many processes which may lead to photon 
absorption and emission in a strong magnetic field, like photon splitting or pair creation via 7 -I- 5 
(resulting then in pair annihilation, and photon energy downgrading via Compton recoil), which 
make photon number conservation unlikely. At the same time, there are important radiation trans- 
fer effects taking place, which obviously tend to favor flatter spectra at low energies; also, all rates 
for photon emission and absorption de part from their values in the absence of magnetic field. De- 



tailed computations ( (lLyubarskyll2002l) ) show that the ensuing spectrum is flatter, not steeper, than 
a blackbody; and we remark here that, despite many calibration uncertainties, fits to the spectrum 
seem to favor spectra flatter than a blackbody, with a preference for thermal bremsstrahlung of 
temperature T = 30 keV. 

We have presented our argument by stressing its independence fr om the estimates of the dipole 
field from p ulsar spin-down, and from Thompson and Duncan's (Thomps on & DuncanI (ll993L 



1995L Il996l) ) model, but it should be obvious that the limit in eq. [13] is perfectly consistent with 



the previous measurements, and the model itself. 
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